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Abstract 

In this paper we discuss the Dirac Oscillator wave equation in terms of a pseudoclas- 
sical language, using Grassmann variables to describe the internal degrees of freedom of 
the oscillator. Regarding the original wave equation as a classical constraint, we use the 
theory of constrained systems, to develop a reparametrization invariant lagrangian, which 
is the pseudoclassical equivalent of the quantum case. The consistency of the Hamiltonian 
formalism and the quantization procedure are also analized. 


1. Introduction 

As is well known, in the second decade of this century Dirac developed the square 
root method, to analyze the internal spin degrees of freedom in quantum mechanics. Dirac 
accomplished his task by means of a Clifford algebra for these degrees. Altought at that 
time the Grassmann algebras were already known, there existed no classical counterpart 
available for his approach. In present day point of view, however, we know that this 
old problem can be formulated using anticommuting, fermionic variables, to reproduce 
the behavior of the spin degrees of freedom at the classical level. Since the Grassmann 
variables have no direct physical meaning, the theories formulated with them are usually 
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called pseudoclassical. 

One important point of the approach formulated above is that, associated with the 
anticommuting variables and involving also the rest of the non-Grassmann dynamical vari- 
ables of the theory, (called bosonic variables), there exist a supersymmetric gauge invari- 
ance in the formulation. One of the reasons for this supersymmetry is the fact that any 
quantum wave equation present in the theory is translated at the classical level as a first 
class constraint. According to Dirac conjecture, all first class constraints generate gauge 
transformations, but since in this case the Grassmann and the bosonic variables are mixed 
up by the gauge transformation they become, in fact, supersymmetric. 

This way of reasoning has been analyzed by several authors [1,2, 3, 4] as a pevious step 
to quantization. The idea is in some sense based in Dirac’s point of view that we should 
first try to fully understand a problem at a classical level, and only then try to quantize it 
[5]. One consequence of this procedure is that we can apply it to systems which we don’t 
fully understand, for instance in the case of two body, or more, relativistic wave equations 
[1]. The interesting point here is that the time evolution of the dynamical variables are just 
the Euler Lagrange equations, which in principle are known, thus solving the dynamics of 
the problem, at least at the classical level. 

2. The Dirac Oscillator 

Let us begin with a simple introduction to the Dirac Oscillator. Some years ago, 
Moshinsky and Szczepaniak introduced a relativistic wave equation linear in momentum 
and in position which has an harmonic oscillator spectrum plus a strong spin-orbit coupling 
term [6]. This equation is obtained by the replacement of the momentum of the particle 
in the Dirac equation! by 


p — > p — imu r/9, 1) 

where p is the momentum, m the mass of the particle and r is its position, u> is the 
frequency of the oscillator, and 0 the Dirac 7 0 matrix. Taking advantage of the frame 
dependence vector u ,L , it is easy to show that the Dirac oscillator equation can be put in 
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the manifestly covariant form 

(6 • P + md 5 )ip = 0, 

where 

P M = p M — 2imu>x^(u • 6)9 5, 

and 

x^ = x M + (u • x)^ 


2 ) 


3. a) 
3.6) 


The operators 6 ** and #5 , are expressed in terms of the Dirac matrices in the following way 


^=^7s 7 m , » = 0,1, 2, 3 

05 =^7s, 


3.c) 


where 75 = i7°7 1 7 2 7 3 . We use natural units in which Ti = c = 1 and our metric is 
given by (?/,,„) = diag.( — 1, 1, 1, 1). We recover the original Dirac Oscillator in the frame in 
which ii # * = (1,0, 0,0). The solution, spectrum degeneracy, hidden supersymetry and other 
important properties are discussed in [8,9], (and references therein). 


3. Pseudoclassical description 

Now, the idea is to reformulate the problem in a pseudoclassical language. This is 
done in a natural way by translating Eq. 2) into a first class constraint 


J = 6 ■ P + m0 5 0, 4) 

where ~ means weak equality and the dynamical variables become pseudoclassical ones 

5. a) 

{$5, ^5} = *• 

Of course we know that 

{x»y} = ti**. 5.6) 

The Poisson braket of the first class constraint J with itself, is the classical equivalent 
of the square of the Dirac Oscillator. In this way, thus, we generate another constraint 
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which could be called the pseudoclassical analogue of the Klein Gordon equation. This 
new constraint is constructed in such a way to be again a first class, altought secondary, 
constraint. Of course, if there is any second class constraint we should replace the Poisson 
braket for the Dirac one, in order to get rid of them. 

In our case the superalgebra obeyed by the constraints is given by 

{ J, J) = iH = i(P 2 + m 2 - AimwX) ~ 0, 

0 , 6 ) 

{H,H} =0. 

In this equation, A is given by 

A = (9 • xj_)u • (mO — 0 5 P). 7) 

Notice that the first of this constraints is precisely the pseudoclassical Klein Gordon equa- 
tion. Is in this sense that we say that we translate the square root method into a classical 
language. We also note from Eq. 7) that H and PL are first class. 

Following the procedure described in Ref. [7], we can construct the Lagrangian of 
the problem. Since the whole dynamics of the theory is contained in the constraints, the 
Hamiltonian of the system is a linear combination of them 

H = NH + iMJ = 0, 8) 

where N and M are gauge fixing parameters. This in turn means that the Hamiltonian is 
weakly zero, implying a reparametrization invariant Lagrangian as a consequence. Accoi- 
ding to Ref. [7] the whole action is given by 

5 = f dr{— my/— z 2 [l — 2iu>(9 ■ ri)(u • 0)] -f + ^5^5] 

Jt2 , , Q) 

- 2 imu)(9 • x±)9 5 {u ■ z) + 2 imuj(z • xx)(u^)^ 5 

— 2 mu>M(9 • xj_)(u • 9)9 5 — imM9$. 

Where z» = - iM9» . 

It is not hard to prove that this action is the correct answer to our problem. From 
the Hamiltonian formalism, we know that the time evolution of any dynamical variable F 
is given by 

F = {F,H). 10) 
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Hence, for the dynamical variables in our case we obtain that 
i" = iM9 M -I- 2 NP 11 + AimuN (8 • r i )^ 5 u'' 

0" = 4m 2 o;iV[x^(u • 9) - (8 ■ xjju' 1 ] + MP" 11) 

8 5 = mM — 2 imMu>(9 ■ x±)(u ■ 8) — 4 mu>N(x± ■ P)(u ■ 8) + Amu>N(8 • x±)(u ■ P ), 
which sure precisely the Euler- Lagrange equations of the action 9) provided 

N(r) = — [1 + 2iu(8 • x±)(u ■ 0)]. 12) 

im 

This result proves the complete agreement between the Hamiltonian and Lagrangian for- 
malisms, and as a result we are confident that action 9) is the corect answer. What Eq. 12) 
tells us is that we can specify the gauge by giving a value to \/—z 2 , usually y—z 2 = —1. In 
the same way we can construct the supergauge transformations for each dynamical variable 
F by means of the eq\iations 

8F= {F,e a (r)<U, a = 1,2 13) 

where <j> a represent our two constraints, and e a (r) are two time dependent infinitesimal 
parameters. The result is too long to be given here, (see reference [7]), the only point we 
want to remark here is that, as we already mentioned, they express the full dynamics of 
the theory, as is suggested by the comparisson of Eqs. 10) and 13). 


4. Conclusions 

We note form Eq. 3. a) that the Dirac Oscillator interaction term is ^-dependent. In 
cases like this, the quantization procedure should be done carefully, since some properties 
of the Grassmann variables changes radically when quantized. For example, the 05 variable 
has the property that 05 2 = 0 at the classical level, but at the quantum level (Eq. 3.c) 
0 5 2 = —1/2. Thus if we consider for instance, the Taylor expansion of a 05 dependent 
function, we obtain different results depending wether the quantization is done before or 
after the series expansion. 

In the case of the system studied here, if we put in Eq. 6) the square of given by 
3. a) and proceed to quantize by means of definitions (3.c), we obtain a wrong result. The 
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central point here, is that we should regard P* 1 as a basic quantity to quantize. If instead 
of developing the square of P M , we first quantize Eq. 6) and regard 3. a) as a quantum 
operator identity, we obtain a complete compatibilty with the Klein Gordon wave equation 
associated with the Dirac Oscillator. Finally, we would like to remark that our approach 
could be useful to problems that are not fully understood at the quantum level, but that 
have 0-dependent interaction terms, such as the aiforementioned two body relativistic wave 
equations and some supergravity theories [7]. 
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